Chirality-spin separation in the Hubbard model on 

^ . the kagome lattice 
(N 

^ . Masafumi Udagawa and Yukitoshi Motome 

' Department of Applied Physics, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, Japan 



C/3 



E-mail: udagawa@ap . t .u-tokyo .ac.jp 



Abstract. Effect of geometrical frustration in strongly-correlated metallic region is studied 
for the Hubbard model on the kagome lattice at half filling by a cluster extension of the 
^ ' dynamical mean-field theory combined with a continuous-time auxiliary-field quantum Monte 

Carlo method. We find that the electron correlation enhances the spin chirality in both vector 
and scalar channels. The chirality grows as decreasing temperature and exhibits a peak at a low 
temperature, indicating a new energy scale under strong correlation. The peak temperature is 
, considerably lower than that for the local spin moment, namely, the characteristic temperatures 

for the chirality and the local moment are well separated. This is a signature of separation 
between spin and chiral degrees of freedom in the correlated metallic regime under geometrical 
. frustration. 

o 

1. Introduction 

J> ' For a long time, considerable attention has been focused on the role of geometrical frustration in 

^ . localized spin systems. It has been discussed that large degeneracy due to geometrical frustration 

\^ I opens a way to realize non-trivial ground states, such as spin-liquid state by suppressing 

Cf^ ' conventional long-range ordering. It has also been clarified that geometrical frustration can 

. induce an exotic long-range order composed of some higher rank objects, such as spin chirality 

O I [3] and spin quadrupole [4j. 

■ On the other hand, much less is known about the effect of geometrical frustration in itinerant 

electron systems. Recently, several interesting phenomena have been found in itinerant systems 
and an importance of the geometrical frustration has been pointed out. A typical example is 
^ . heavy- fermion behavior found in transition metal compounds such as LiV204 [5], Y(Sc)Mn2 

^ I [6] and /3-Mn [7j. In contrast to the heavy-fermion rare-earth compounds, these materials 

have no explicit localized moments, and therefore, it is difficult to explain their heavy-fermion 
behavior by the conventional Kondo mechanism based on the interaction between conduction 
electrons and localized spins. A candidate for the origin of the heavy-fermion behavior is a 
common feature among these transition metal compounds — geometrical frustration in the 
underlying lattice structure. In fact, in the related insulating materials, it has been reported 
that magnetic transition temperatures are strongly suppressed compared to the Curie- Weiss 
temperatures because of the frustration [8]. Furthermore, characteristic magnetic fluctuations 
under a suppression of long-range ordering are observed in these materials, and their relation to 
the heavy-fermion behavior has attracted much interests [9l \T0\ [TT] . 

In general, in a weakly-correlated metal, the electronic state is well described by a Slater 
determinant of single-particle states which extend over the entire system. In this case, physical 



properties of the system will be rather insensitive to the local lattice structure, i.e., whether or 
not the lattice is geometrically frustrated. The issue is how this picture is modified when the 
system enters into a strongly-correlated regime where the electrons tend to be localized due to 
strong electron interaction. In particular, in the vicinity of the Mott transition, spin moments 
grow with antiferromagnetic correlations between neighboring sites, and hence, it is expected 
that these moments suffer from the frustration and lead to some exotic behavior as seen in the 
localized spin systems. It is intriguing to examine how the geometrical frustration affects the 
electronic state through the interplay between spin and charge degrees of freedom. 

Motivated by these considerations, in this paper, we investigate the effect of geometrical 
frustration in correlated metallic region for a simple model, the Hubbard model on the kagome 
lattice. In particular, we will focus on the behavior of spin chirality in comparison with that of 
the original spin degree of freedom. 

2. Model and Methods 

As a minimal model including both geometrical frustration and electron correlation, we consider 
the Hubbard model on the kagome lattice, whose Hamiltonian is given by 



in the standard notations. We consider only the nearest-neighbor site hopping on the kagome 
lattice shown in the left panel of Fig. [H and set t = 1 as an energy unit hereafter. The chemical 
potential fi is controlled so that the system is at half filling (one electron per site on average). 
The Mott transition and magnetic fluctuations near the transition have been recently studied 
for this model by several theoretical methods \12\ [T3\ [T3] . 

To investigate the correlated metallic region of the model (1), we use a cluster extension of the 
dynamical mean-field theory jT5|, which enables to describe the Mott transition with including 
spatial fluctuations within a finite-size cluster. In this method, the original lattice problem is 
mapped to a cluster impurity problem: In the present study, in order to take account of the effect 
of geometrical frustration, we consider a 3-site cluster (Fig. [1]). We extended the cluster up to 
12 sites and confirmed that the spin correlation hardly develops beyond the nearest-neighbors in 
the parameter range of the following calculations (not shown). To solve the impurity problem, 
we employ the continuous-time auxiliary-field quantum Monte Carlo method [16j, which requires 
much smaller computational cost compared with the conventional Hirsch-Fye algorithm. 

In the next section, we will show the results for the density of states, the spin chirality and 
the local spin moment. The density of states (DOS) p{uj) is calculated from the one-particle 
local Green's function, Gjair) = — (TrCj(j{T)c'j^{Q)) , (j = 1 — 3 corresponds to the cluster sites in 

Fig.[T]) by using p{oj) = —^lm.Gja{oJ+i6). We use the maximum entropy method for the analytic 
continuation. Note that p{uj) does not depend on either j or a within the paramagnetic solution 
assumed here. For the calculations of magnetic properties, we define the local spin at cluster 
site j as = i Yls s' (^\si^^)ss'Cjs' ■, with the v-ih component of Pauli matrix, {u = x, y, z). 
By using Sj, we calculate the vector spin chirality, = ^^(^i x S2 + S2 x S3 -|- S3 x si), and 

the scalar spin chirality, = (si x S2) • S3. In the following, we present the squared moments 
of the vector chirality K"^ = (K^) and of the scalar chirality = (K^) as well as the local spin 
moment = (sj) (s^ is independent of j). Here, the bracket represents the statistical average 
taken in the grand-canonical ensemble at temperature T. s^, and are two-, four- and 
six-body equal-time correlation functions, respectively, and obtained by applying the Wick's 
theorem for each Monte Carlo sampling. 
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3. Results 

First, we present the result of DOS at [/ = 6 and T = 0.05 in Fig. [2j For comparison, we also 
show the non-interacting DOS at C/ = 0. At C/ = 0, p{oj) structures, two van-Hove singularities 
(a; ~ —0.47 and —2.47) and a flat band (cj ~ 1.53). At f7 = 6, the system is still in the metallic 
region with finite DOS at the Fermi level w = 0, whereas p{uj) shows several characteristic 
features of strong correlation effect. One is the strong renormalization of the energy scales; the 
three peaks around cj = are considered as renormalized structures of the van-Hove singularities 
and the flat band. Another feature is the formation of the Hubbard bands which are observed 
as two broad humps at w ~ ±3. We note that the result is consistent with the previous one 
obtained by a similar method but with using the conventional Hirsch-Fye algorithm [14J. The 
previous study revealed that the critical value of U for the Mott transition is about 8.2. Below 
we show the results for the spin chirality degrees of freedom in this correlated metallic region at 
U = 6. 





Figure 1. Schematic picture of the 
mapping of the Hubbard model on the 
kagome lattice to an effective cluster- 
impurity Anderson model. 



Figure 2. DOS at [/ = 6 and T = 0.05 
(red) in comparison with the result at 
U = (black). We set the chemical 
potential to be cj = 0. 



Figure El shows the temperature dependence of the spin chiralities and K'^. At [7 = 0, 
both and are small and featureless, in spite of the peculiar u dependence of the DOS. 
The deviation from the high temperature limits [K^{T — > oo) = 1/24 and K'^iT — )• oo) = 3/256] 
is within 0.03 % for the entire temperature range. In contrast, at [/ = 6, both i^^ and 
are largely enhanced compared with the non-interacting case. Furthermore, as decreasing T, 
they grow gradually and exhibit a peak at a low temperature ~ 0.3. This temperature 
characterizes a new energy scale related with the spin chirality degrees of freedom. 



0.15 



0.05 





' k/(U=6) X 




K5^(U=6) □ 


i **■ . 


k/(U=0) 


m 

m 


K32(U=0) 




m 




'# 




m 


• " ■ ■ - . 













0.0 0.4 



0.8 1.2 1.6 

T 



1.75 



1.25 



1.00 




0.0 0.4 



1.06 



1.00 



0.97 



Figure 3. Temperature dependence of 
Kl and K'i &t U = Q &nd U = 0. 



Figure 4. Temperature dependence of 
Kl, K"^ and at [7 = 6. Ah the data 
are normalized by the value at T = 2.0. 



In Fig. m we compare this behavior of the chirahties with the local spin moment s^. also 
shows a broad peak, but the peak temperature T* ~ 0.9 is considerably higher than T^. In the 
intermediate temperature < T < T* , the chiralities are enhanced while the spin moment 
is suppressed as decreasing T. This contrastive behavior is surprising because of the following 
reasons. First of all, by definition, the spin chiralities are composed of the products of spin 
operators. In fact, their behavior is tightly related with that of the spin moment in most of 
the localized spin systems. Furthermore, in the present itinerant system, the decrease of is 
anticipated to reduce the spin chirality moments since it corresponds to the increase of doubly- 
occupied or empty sites: The spin chiralities have largest values when they operate on the 3-site 
triangle which includes no doubly-occupied or empty site. Therefore, the apparent separation 
between the spin and chirality found in Fig. H] is highly non-trivial. 

4. Summary and Concluding Remarks 

We have studied the spin and chirality moments for the Hubbard model on the kagome lattice at 
half filling by applying the cluster dynamical mean-field theory combined with the continuous- 
time quantum Monte Carlo method. In metallic region under strong electron correlation, we 
obtained the following results. (1) The spin chirality is strongly enhanced in both vector and 
scalar channels, compared with the non-interacting case. It grows as decreasing T and exhibits a 
peak at a low temperature T^, implying a new energy scale due to the electron correlation. (2) 
The local spin moment also shows a peak, but the peak temperature T* is substantially higher 
than T^. In the temperature range between and T* , the spin and chirality show opposite 
temperature dependences. These results suggest a separation of the spin and chirality degrees 
of freedom in the correlated metallic region in this strongly frustrated system. 

The origin of the chirality-spin separation deserves further study. It is also interesting to 
clarify how the electronic state is affected by the appearance of the new energy scale related 
with the chirality. Another issue is a possibility of some instabilities in the chiral sector by 
allowing symmetry breaking within the cluster dynamical mean-field framework. Systematic 
study in a wide region of U and T is now in progress. 
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